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Abstract 

Lensing effects on light rays from point light sources, such like Type la super- 
novae, are simulated in a clumpy universe model. In our universe model, it is assumed 
that all matter in the universe takes the form of randomly distributed objects each 
of which has finite size and is transparent for light rays. Monte-Carlo simulations are 
performed for several lens models, and we compute probability distribution functions 
of magnification. In the case of the lens models that have a smooth density profile or 
the same degree of density concentration as the spherical NFW (Navarro-Frenk- White) 
lens model at the center, the so-called gamma distributions fit well the magnification 
probability distribution functions if the size of lenses is sufficiently larger than the 
Einstein radius. In contrast, the gamma distributions do not fit the magnification 
probability distribution functions in the case of the SIS (Singular Isothermal Sphere) 
lens model. We find, by using the power law cusp model, that the magnification prob- 
ability distribution function is fitted well using the gamma distribution only when the 
slope of the central density profile is not very steep. These results suggest that we may 
obtain information about the slope of the central density profiles of dark matter halo 
from the lensing effect of Type la supernovae. 
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§1. Introduction 

Type la supernovae (SNe) are very useful tools to investigate our universe. For exam- 
ple, the distance-redshift relation obtained as a result of the observations of these strongly 
suggests the present acceleration of the cosmic volume expansion of our universe. 1 ^ 4 - ) To 
confirm this indication, further projects, such as the ESSENCE project at NOAO0 the 

; are now 



Large Synoptic Survey Telescope**- 1 ! and Supernova/Acceleration Probe (SNAP) 
active or planned. 

Current observational data points on the distance-redshift plane are somewhat scattered, 
and this is an origin of the error in estimating the cosmological parameters. One of the 
reasons for this distance dispersion is gravitational lensing effects due to the mass inhomo- 
geneities in our universe. Although the dispersion due to gravitational lensing is relatively 
smaller than other effects at low redshift, it may become more prominent and comparable to 
the intrinsic dispersion at high redshift z > I.2. 5 )' 6 -' In this sense, the mass inhomogeneities 
are obstacles to the determination of the cosmological parameters through the observation 
of Type la SNe. On the other hand, the mass inhomogeneities contain rich information 
about the physical process of the evolution of our universe. Thus, the mass inhomogeneities 
themselves are very significant subjects in cosmology. 

We can understand the property of the mass inhomogeneities by comparing the observa- 
tional data about type la SNe with theoretical predictions obtained through the investiga- 
tions of the gravitational lensing effects on the light rays from them. This is just the purpose 
of this paper. 

There are many works about gravitational lensing related to SNe: effects on SNe obser- 
vations, 6 )~ 15 ) their availability for investigation of our uni verse, 16 )~ 22 ) the extraction of the 
evidence for lensing effects from observational data 23 - 1-26 '' and others. 27 -*' 28 )> 30 )~ 36 ) 

In this paper, we particularly focus on the effects of small-scale structures. The smaller 
scale inhomogeneity may affect the apparent magnitude of sources owing to gravitational 
focusing. The mass scale of the inhomogeneities, which is considered in this paper, is between 
10~ 2 M Q and 1O 1O M . There are several works related to this subject. 16 M 8 )> 27 )> 28 ) i n these 
works, compact lens objects were mainly studied, whereas, in this paper, we consider the 
extended lens objects. Even in the case of the extended lens objects, the strong lensing 
effects that cause multiple images are of considerable importance if the density profiles of 
the objects are significantly steep at the center. 



*) http:/ /www. ctio.noao.edu/wproject/ 
**) http:/ /www. lsst.org/lsst_home.shtml 
***) http://snap.lbl.gov 
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We consider the flat ylCDM model with Qao — 0.71, i? m o = 0.29 as the starting point. 
In the cold dark matter cosmology, the first objects to form are of subgalactic size. Then, 
larger structures form through tidal interaction and mergers of smaller objects. Therefore, 
the dark matter may form numerous clumps. For simplicity, we assume that all matter in 
the universe takes the form of randomly distributed objects, each of which has finite size and 
is transparent for light rays. We investigate the magnification effect due to the gravitational 
lensing on light rays from the point sources such as Type la supernovae. The light rays from 
distant supernova events suffer multiple gravitational lensing effects from the clumps in the 
universe. The propagation of the light rays is thus stochastic owing to the randomness of the 
distribution of the clumps. In this paper, we focus on magnification probability distribution 
function (MPDF). We introduce a simple method to simulate multiple gravitational lensing 
effects in the clumpy universe model, which is based on our previous study. 38 -* We perform 
Monte-Carlo simulations for several lens models, and we investigate the dependence of the 
MPDF on lens models of dark matter clumps. 

Throughout this paper, it is assumed that light sources are pointwise. This assumption 
is valid only when the linear extent of a source is much smaller than the Einstein radius 
of each lens object. The linear extent of the SNe is of the order of ~ 10 15 cm, which is 
the typical radius of SNe at the peak brightness. The Einstein radius of a lens object of 
mass Ml is of the order a/ GMlD/c 2 , where D is the distance to supernovae, G and c are 
the gravitational constant and speed of light, respectively. Therefore, we assume that Ml 
satisfies \JGM L D / c 2 > 10 15 cm. Because in the case of high-redshift supernova events, 
D ~ lGpc, this assumption leads 

M L > 1(T 2 M . (1-1) 

Our method is complementary to the methods using the N-body simulation. In the 
N-body simulation, the nonlinear time evolution and the spatial correlation of the dark 
matter distribution can be studied. Very large scale N-body simulations are possible now. 
Nevertheless, the questions that can be definitively answered with N-body simulations are 
still limited by finite resolution in mass and in distance. 29 ^ Since the small-scale structure 
of the dark matter is essential for the gravitational lensing effects, the simulations using 
analytic density profiles are important to confirm the results from the N-body simulation, 
and to investigate new effects that are not found in the N-body simulations. 

This paper is organized as follows. In §2.11 we briefly review the gravitational lensing. 
§2.21 is devoted to describe the clumpy universe model. In §2.31 we show a definition of 
magnification and the relation between the distance and the magnification. In £j3j numerical 
calculation and our results are illustrated. Finally, §U is devoted to the conclusions and 
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summary. Throughout this paper, we use the unit of G = c = 1. 



§2. Gravitational lensing in a clumpy universe model 



2.1. Lens equation 

Before describing our settings and calculation method, let us briefly review gravitational 
lensing. In this paper, we focus on the situation in which thin lens approximation 37 ) is valid. 
The thin lens models are characterized by the surface mass density 2?(£), where £ is the 
impact vector (see Fig. [I]). By using the surface mass density the bending angle 

vector is given as 37 -* 



For convenience, we consider the "straight" line A from a source to the observer and define 



Observer • * * — 

Fig. 1. Geometry of gravitational lensing due to a point-mass lens. The impact vector £ represents 
the relative position of a light ray on the lens plane to the lens position, and a is a vector whose 
norm is equal to the bending angle. 

the intersection point of A with the lens plane as the origin of the lens position £ and the 
ray position 7. In the geometrically thin lens approximation, we can regard the vectors £, 
£, and 7 on the lens plane. 37 ) 

We define Ds, Di, and Dis as the angular diameter distances from the observer to the 
source, from the observer to the lens, and from the lens to the source, respectively. Using 
simple trigonometry (see Fig. [TJ, we find the relation between the source position 77 and £, 




(2-1) 




77 



(2-2) 
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From Eq. ( \2-2L we have 

= ^7 " Dlsol^ - (2-3) 

by using 7 instead of £ and 77. 

2.2. Clumpy universe model 

The clumpy universe model in which we calculate the lensing effects is basically the same 
as that in Ref. 38). In this clumpy universe model, it is assumed that all matter takes 
the form of randomly distributed objects and its global property is well described using the 
Friedmann-Lemaitre (FL) universe whose metric is given as 

ds i = + a 2( t ) (-^j^i + , (24) 

where K — 1, 0, and —1, and dQ 2 is the round metric. Hereafter, we will refer to this FL 
universe as "background universe" . 

We assume that the comoving number density p n of the lenses in the clumpy universe is 
given as 

_ a 3 p _ 3f2 m0 HQ 
Pn ~ M L ~ 8vrM L ' 

where p, i7 m0 , and H are the average mass density, the present values of the total density 
parameter, and the Hubble parameter, respectively. We consider a past light cone of the 
observer at r = in the background universe, which is parametrized by the redshift z of its 
null geodesic generator. The comoving volume AV of a spherical shell (see Fig. [2]) bounded 
by r = r(z) and r(z + Az) on this light cone is given as 

47rr 2 dr . . 
AV = —Az. 2-6 

VI + Kr 2 dz 

Therefore, the number of lenses in this shell is given as 

AN = Pn AV= 3 -^SK d ^Az. (2-7) 
2M L Vl + Kr 2 dz 

We assume dr/dz is the same as the background value, 



dr J_ / l + H 2 f2 K0 (l + z) 2 D 2 FL (z) 
dz H V a m0 (l + zf - Q KQ {1 + z) 2 + O A0 



(2-8) 



where i?Ao an d Dp^z) are the normalized cosmo logical constant and the angular diame- 
ter distance from the source of the redshift z to the observer in the background universe, 
respectively, and ^?ko = ^mo + ^ao — 1- 



5 



Let us consider a lens plane at the redshift z. We define y as 

c 



y 



where £o is the Einstein radius given as 



DrD 



(2-9) 



(2-10) 



y represents the lens position normalized with the Einstein radius. In addition, we write 
the absolute values of y and £ as y and (, respectively. The average number of lenses in the 
region [y, y + Ay) within the domain [z, z + Az)(see Fig. |2J) is given as 



47ra 2 r 2 



-AN, 



(2-11) 



where we have assumed ( <C r. This value is equivalent to the probability p(y,z)Ay for a 

★ 

Source 



Line A ! 




Line A 



r(z + Az) --- 




r{z + Az) 



Fig. 2. Figures for E fO 

light ray to receive lensing effects with impact parameter [y, y+Ay) in the domain [z, z+Az). 
Substituting Eqs. §FQ) and (JMD into Eq. (j2TT|) and using Eq. fT2~TUD . we have 

3H ^mo(l + zf 



p(y,z)Ay 



DrD 



L^LS 



yAyAz. 



(2-12) 



v/W + ^) 3 - ^o(l + + Rio ^5 

It is seen from Eq. (12-121) that this probability is specified using the model of the background 
universe. 

2.3. Magnification and distance 

In this subsection, we discuss the definition of magnification due to lensing and the 
definition of distance in the clumpy universe model. 
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First, we consider a point light source that is observed as a single image in a clumpy 
universe. Magnification ji for a source is defined as 

S 

»=-=-, (2-13) 
>->o 

where S is the flux actually observed, whereas So is the fictitious flux that would be emitted 
by the source with the same luminosity and redshift but without lensing effects. If the light 
ray bundle passes through a spatial section with an area AA orthogonal to the ray bundle, 
using the conservation law of photon number, 37 ) we find 

S AA 

"=sr:zr (2 ' 14) 

where AA is a sectional area orthogonal to the ray direction in the fictitious propagation 
process. 

The corrected luminosity distance is defined as 

where AQ S is the solid angle that subtends a light ray bundle from the point source. The 
relation between the angular diameter distance D and the corrected luminosity distance D' 
is given as D' — (1 + zs)D, where z s is the source redshift. 37 ) Let D denote the angular 
diameter distance between the source and observer in the case without any lensing effects. 
Since 



An* 



= (l + z s )D , (2-16) 



we obtain 



0' = f^(l+^)A, = (2-17) 

\AQ S J 

Therefore, we can define the observed angular diameter distance Da as 

D A = ^. (2-18) 
v A* 

Next, let us consider the case of multiple images. The magnification /i^ of the p-th 
image of a source is defined as 

m = (2-i9) 
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where S( p ) is the flux of the p-th image. In this study, we assume that multiple images can 
not be distinguished from each other owing to the limitation of the resolution, and we can 
observe only the total flux of light rays. Then the observed magnification p is given as 

^ = ^ = E^)» ( 2 - 2 °) 

>->o 

p 

and we define the observed angular diameter distance as 

D A = D ° . (2-21) 

y E P Hp) 

In terms of the angular diameter distance, the magnification is written as 

D 2 

= (2-22) 

A 

The angular diameter distance in gravitational lensing is discussed for many years. 45 -*' 46 ) 
The angular diameter distance in the FL universe, Dfl, and the Dyer-Roeder distance, 40 -*' 41 -* 
Ddr, are commonly used. The Dyer-Roeder distance is the distance in the under dense region 
in the inhomogeneous matter distribution. It depends on the smoothness parameter a, which 
is the ratio of the smoothly distributed matter density except for clumps to the mean energy 
density p for all matter in the universe. When a — 1, the universe becomes the FL universe 
and Ddr agrees with Dfl- When a — 0, all matter takes the form of clumps, and in such a 
case, the Dyer-Roeder distance is the distance in an empty region in the clumpy universe. 

In our work, unlensed angular diameter distance Dq should be appropriately chosen so as 
to be consistent with the situations of our present interest as below. 51 -* We are investigating 
the lensing effects of SNe, whose linear extent is much smaller than the scales in which 
density inhomogeneities are linear or quasi-linear. Then, the matter distribution might be 
highly clumpy, and the under dense region might be almost empty. 

In this situation, we expect that the focusing effects the light ray receives are minimal, 
and that the angular diameter distance agrees with the Dyer-Roeder distance, because the 
light rays do not receive any lensing effects between the lens planes. This is consistent 
with the numerical simulations by Holz & Wald 27 ) and Kozaki. 52 ) Holz & Wald 27 -* discussed 
the angular diameter distance in the universe that consists of many small balls. Kozaki 52 ) 
calculated the distance-redshift relations in the Swiss-cheese universe model, that is an exact 
solution of the Einstein equation. In this paper, we adopt the Dyer-Roeder distance D DR 
with a = as D and for the distance between each lens plane which has positive definite 
surface mass density Hereafter, we drop "with a = 0", and simply call this distance 

the "Dyer-Roeder distance Ddr". 
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In Refs.49), 50), it is shown that the observed distance after the multiple lensing effects 
sometimes becomes longer than the DR distance. In such cases, a becomes negative. In our 
simulation, we do not treat such a case and use the Dyer-Roeder distance as the angular 
diameter distance. To understand such situations, we need further studies. 

If inhomogeneities are so small that these are regarded as linear or quasi-linear per- 
turbations in the homogeneous and isotropic universe, perturbative treatments seem to be 
better in gravitational lensing. In such perturbative treatments, the unlensed distance D Q 
is the angular diameter distance in the FL universe Dfl- Light rays are perturbed by the 
lensing effects of the inhomogeneity. In this sense, the (de)magnification is often defined as 
fi := D 2 FL /D 2 A in weak lensing analyses and ray shooting calculations. 42 )~ 44 ) In their works, 
the surface mass density of each lens plane is not positive definite, namely, the under dense 
regions have negative mass density. The relation between /i and pL is given as 



3.1. Numerical methods and the density profiles of lens models 

In this study, we use the numerical method proposed by Rauch. 7 ) We compute lens 
mappings using the multiple lens plane method. ^ We ignore the lensing effects coming 
from clumps far from the ray. We search for the position of images on the field of view using 
the Newton-Raphson method. Summing up the magnifications of all images, we obtain the 
total magnification. The concrete description about the method is in Ref.38) and we review 
it in Appendix lAl 

In the simulations, we assume that clumps in the universe are transparent to light rays. 
We also assume that clumps are axisymmetric, and the axes of symmetry of clumps are 
identical with the line of sight. Each of these clumps has the mass Ml and the "size" R. 
We consider the following three lens models: 

(a) Homogeneous disk: 



FL 



(2-23) 



DR 



§3. Numerical calculation and results 




for f < R, 
for f > R. 



(3-1) 



*) Numerical calculations have been performed also for homogeneous sphere:I7(£) = §§jp (l — £ 2 /R 2 ) 1 ^ 2 
and power low tail model:Z'(^) = ^2^^2/^.2)2 ■ The results are qualitatively the same as those of the lens 
models (a) and (b). We do not show these results in this paper. 
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(b) Log cusp model: 

27(0 = i ^ ~ ? "~ * ^ (3-2) 




(c) l/£ cusp model (Singular Isothermal Sphere): 

f 9%T for £ < 
^(0 = 2 ^ n ? , ' (3-3) 
I tor 4 > ft. 

The degree of the density concentration at the center of the log cusp model (b) is equiv- 
alent to the spherical Navarro- Frenk- White (NFW) 39 ^ lens model in which the mass density, 
p, is inversely proportional to the radius. 

The surface mass density of all our lens models can be written in the form, 

where F(R/£) is a function of R/£, £ is the distance from the axis of symmetry on the lens 
plane, and 27(0 * s the surface mass density. We assume that the total mass of each lens 
object Ml = 27r J" °° £27(£)^£ is finite. If we use Eq. (13-41) and use the calculation method 
written in Appendix [Aj we find that the magnification distribution does not change if we 
scale the parameter in the manner R — > CR and Ml — ► C 2 Ml, where C is an arbitrary 



constant. This means that the magnification depends only on the parameter Rj \JMl- A 
proof of this fact is given in Appendix [A] 

3.2. MPDFs 



We first consider the cases in which Rj \JMl =constant, and investigate the dependence 
of MPDFs on lens models and on the value of Rj \/M£. The magnification is stochastic 
owing to the randomness of the lens distribution. We numerically generate 100,000 samples 



of magnification for each lens model characterized by the parameter Rj \JMl- In our cal- 
culations, we use the same realization of the lens distribution for each model. The MPDFs 
given as results of these calculations are shown in Figs. [3H5J 

In Fig. [3j MPDFs for case (a) are depicted. There are a few frames in which two or three 
peaks appear in an MPDF. The largest peak and the second or third largest peak in the 
case of (zs, R\J H /M L ) = (0.4, 6), (1.2, 2), and (2, 2) can be explained as the effect of the 
small lensing probability. When the lensing probability is relatively small, most of the light 
rays do not experience the lensing. In this case, it is thus expected that the MPDF has a 
peak at 1. This is because, in our definition of the magnification, when a light ray does not 
experience the lensing, the magnification is 1. When the light rays experience the lensing 
effect only once, the magnification typically becomes p ~ R 4 /(R 2 — MlDs) 2 as shown in 
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Appendix iBl As a result, we find the second peak around fi ~ R A /(R 2 — MlDs) 2 , which 
can be seen in Fig. [3j 

In Figs. [3]l5l we find that, in many cases, MPDFs can be fitted using the gamma distri- 
bution /(// — l;k,9) defined as 



p-tn-i)/« 



f^-l- k ,9) = {»-lY^- §j —^, (3-5) 

where k > and 9 > are called the shape parameter and scale parameter, respectively, and 
r(k) is the gamma function. In these figures, we depict the gamma distributions that fit 
the MPDFs. The parameters of the gamma distributions are determined by the maximum 
likelihood method as follows. We assume that the numerically generated magnifications are 
independent of each other and identically distributed random variables that obey the gamma 
distribution. The logarithm of the likelihood function £ for A" samples of the magnification, 
Gui, . . . ,/xjv) is given as 

N 

\n£(k,9) = Y,^f(fi t -l,k,9). (3-6) 

i=l 

The shape parameter and the scale parameter are related to the mean /x m and variance 
a of the gamma distribution as follows: 

fi m = k9 + l, (3-7) 
a 2 = k9 2 . (3-8) 

The maximum likelihood estimates of the shape and scale parameters, k and 9, are obtained 
by solving the equations, 



d 



09 







— \nC(k,9) . . = 0, — ln£(fc,0) - . = 0. (3-9) 



k=k,e=e dk 



k=k,e=e 



The means (/i) and variances of the numerically generated MPDFs are plotted as func- 
tions of zs in Figs. [6] and [3 The shape parameter k and the scale parameter 9 of the gamma 
distributions that fit the MPDFs are also plotted as functions of zs in Figs. [8] and [9j As 
explicitly shown in Fig. 0, (//) is almost equal to D 2 DR /D 2 FL . The deviations are smaller 
than 1% in all cases. 

In the case of z$ = 1.2 and z$ = 2.0 for the lens models (a) (see Fig. [3]), the gamma 
distributions fit well the MPDFs for Rj ' s[Ml = 6H~ 1/2 , 10H~ 1/2 but do not for R/^fMl = 

1 /2 

2H L/ . In the case of model (b) and z$ > 1.2 (Fig. d]), the gamma distribution fits 
well the MPDFs for R/VMl > 2H ' . In contrast, the gamma distributions do not fit 
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MPDF at z s =0.4 for the model (a) 



MPDF al z s =1.2 for the model (a) 



MPDF at z s =2.0 for the model (a) 
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Fig. 3. MPDFs for the lens model (a) at z$ = 0.4, 1.2, and 2.0 are shown on the left, center, and 
right lines, respectively. The smooth lines are the gamma distributions that fit the MPDFs. 



well the MPDFs in all cases for lens model (c)(see Fig. [5]). The lens model (c) has the 
steepest density profile at the center among the lens models considered here. The reason for 
the deviation from the gamma distributions seems to come from this steep density profile. 
This is consistent with the fact thatjin the case of the point mass lens model, the gamma 
distributions do not fit the MPDFs. [j This also suggests that, even in models (a) and (b), 
if we set Rj\JM L sufficiently smaller than R/\/M L = 2H 1 ^ 2 , the gamma distributions do 
not fit the MPDFs irrespective of zs- 

1/2 

Let us discuss the z$ dependence of MPDFs. We consider the case of R/^/Ml = 6H 
for the lens model (a). While the gamma distribution fits well the MPDF at z$ = 2, it 
does not at zs = 0.4. The most significant difference between these two cases is the number 
of lensing effects the light rays experience. Hence, our results suggest that, in order for a 
MPDF to be fitted well by the gamma distribution, the light rays need to pass through near 
the clumps frequently. However, to clarify the condition for the realization of the gamma 



*) In the case of point mass lens, the magnification probability behaves ~ /i 3 when [i ^> 1 (see, e.g. 
Ref.53) or 54)). 
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MPDF at z s =0.4 for the model (b) 



MPDF at z s =1.2 for the model (b) 



MPDF at z s =2.0 for the model (b) 
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Fig. 4. 
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Fig. 4. Same as Fig. EJ but for the lens model (b). 



distribution more precisely, much more detailed investigations are required. 

The above results suggest that if lens models have sufficiently gradual density profiles, the 
MPDF is universally well fitted using the gamma distribution. To obtain further evidence 
for this hypothesis, we calculate an MPDF in the case when the value of Rj \/Ml is not 
constant. We assume the lens model (b) in which R^J H Q / of each lens randomly takes 
a value within 6 < R\J H /Ml < 10, and that the redshift of the light source is given as 
z s = 1.2. The MPDF is shown in Fig. EH One can find that MPDF is fitted well by 
the gamma distribution even if the value of R/y/M£ is distributed. This fact might not 
be surprising, because the effect of changing the value of Rj y/Mi, is similar to the effect of 
changing the impact parameter of a light ray from a lens. Both effects produce the change 
in the magnification. The shape parameter and the scale parameter become 9 = 0.041 
and k = 4.5, respectively. These are intermediate values between those in the cases of 
Rj^fMl = QH 1/2 and Rj ^W L = IOHq 12 for the lens model (b). 
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Fig. 7. The variances are depicted as functions of the source redshift. 
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Fig. 8. The shape parameters k are depicted as functions of the source redshift. 



3.3. x 2 -test 

The discussions in §3.21 suggest that we may obtain the information about the slope of 
dark matter density profiles at the central cusp from the Type la supernovae observation. 
To clarify this possibility, we investigate the relation between the goodness of fitting of the 
MPDFs to the gamma distributions and the steepness of the central density profile using 



15 




,. ,,, 1/2 -1/2 
R/M L =6H 




1.6 1.8 2 



V2- jnzj -m 



0.055 

0.05 
0.045 

0.04 
0.035 
CD 0.03 
0.025 

0.02 
0.015 

0.01 
0.005 



R/M L =10H 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1 



Fig. 9. The scale parameter is depicted as a function of the source redshift. 
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Fig. 10. MPDF for the lens model (b) at z$ = 1.2 is shown. Values of Hq/Ml of clumps 
are distributed uniformly within 6 < Ryf Hq/Ml < 10. The smooth lines are the gamma 
distributions that fit the MPDFs. 



the power law cusp lens model denned as 

(2-n)M L 



for £<R, 



= { 2 * R2 \SJ — * - ( 3 .10) 
for f > R. 

This lens model includes the lens models (a) and (c) in §3.21 as the special cases. If we set 
n = or 1 in the Eq. (13-101) . this model becomes model (a) or (c), respectively. 
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We use the value of x 2 as an indicator of the goodness of fitting to the gamma distribution. 
Here, x 2 is defined as follows. The values of the magnification /x of each sample are labelled 
with i (i — 1, . . . , N) so that /ij satisfies the relation, jjli < /!% < . . . < /xjv. We divide the 
region jj, > 1 into p bins. In this paper, we set p = 15. The boundaries of the bins are 
expressed as Vj (j = 0, . . . ,p — 1), which satisfy the relation, is < V\ < . . . < v p _i, where 
uo = 1. We define V\ = (fi m + fi m+ i)/2 so that there are m samples in the 1st bin. In the 
same way, we define v v -\ = (/ijv-m + A*iv-m+i)/2 so that there are m samples in the pth bin. 
The values of Vi from % = 2 to p — 2 are determined by even spacing. Namely, they are given 
as 



Vj=vi + {3 -l)Av, (j = 2,...,p- 1) 
Av = {v p „ 1 -v 1 )/{p-2). 



(3-11) 
(3-12) 



The value of m is determined such that each bin has a sufficiently large number of samples. 

For j-th bin, we have the number of samples in this bin, rrij, and the expected number 
of samples, ej, derived from the gamma distributions. We then define the value of x 2 as 

15 



X 



-2 ^ m i lo S 



8=1 



_ej_ 



(3-13) 



In Fig. [TT], we depict the value of x 2 as the function of the steepness n of the density profile, 
where N = 100, 000, and we have set R/yfMl = WH~ 1/2 , z s = 1.2, and m = 100. The 




Fig. 11. The value of x 1S depicted as the function of power n, where we set R/y/Mi, = 10i^ 
and zs = 1.2. 



1/2 



value of x 2 increases rapidly with n if n > 0.6. 

SNAP will be able to find the order of 1000 SNe per year within the redshift ~ 1. To 
simulate such situation, we divide the 100,000 samples into 100 datasets of 1000 samples, and 
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perform the x 2_ t es t for each dataset. In this case, we set m = 10. When the value of \ 2 is 
larger than a certain threshold, we judge that the MPDF cannot be fitted using the gamma 
distribution. The threshold is determined such that there is only a 5 % probability that the 
dataset generated from the gamma distribution is rejected. The relation between the rate of 
rejection and the value of n is shown in Tabled We find that if n < 0.7, the rejection rate 
is around 10% or less than that. However, if n > 0.8, the rejection rate increases rapidly. 
This result suggests that we may constrain the slope of the dark matter density profile at 
the central cusp from the goodness of the fitting of the gamma distribution to the MPDF. 



Table I. The rate of rejection through the x 2_ test with 1000 samples is shown in this table. 



Power n 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


Rate of rejection 


0.1 


0.06 


0.06 


0.05 


0.12 


0.14 


0.07 


0.11 


0.27 


0.82 


0.97 



§4. Summary and conclusions 

We have studied lensing effects on observations of point sources such as Type la super- 
novae by the Monte-Carlo simulation. It has been assumed that all matter in the universe 
takes the form of randomly distributed objects each of which has finite size and transparent 
to light rays. In addition, we have assumed that each of the lens objects is axially symmet- 
ric along the line of sight. We have found analytically that the magnification probability 
distribution functions (MPDFs) depend on the mass Ml and the size R of a lens object 
only through the form Rf \JMl- We have calculated MPDFs for various values of the source 
redshift and R/\/Ml- We found that the resultant MPDFs can be categorized into two 
groups according to whether the gamma distributions fit well the MPDFs or not. In the 
case of the lens models that have smooth density profiles, the gamma distributions fit well 
the MPDFs if R^J Ho /Ml is sufficiently large. Furthermore, we have found that this result 
holds in the lens model that has the same degree of density concentration at the center as 
the Navarro- Frenk- White (NFW) lens model. In contrast, the gamma distributions do not 
fit well the MPDFs for any value of Rj \JMl in the case of the singular isothermal sphere 
(SIS) lens model. These results suggest that we might be able to distinguish the NFW lens 
model from the SIS lens model on the basis of MPDFs. 

We have shown that MPDF is well fitted using the gamma distribution even if the value 
of R/^Ml is uniformly distributed within 6 < R/^/Ml < 10 in the case of lens model (b). 
This result suggests that MPDF is universally well fitted using the gamma distribution if 
lens models have sufficiently gradual density profiles. To obtain further evidence for this 
hypothesis, we need to clarify the effect of the distribution of R/^/Ml in more detail. This 
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is one of our future works. 

We introduced the power law cusp lens model, and we investigated the dependence of 
the goodness of the fitting with the gamma distribution on the power n, which represents 
the steepness of the cusp at the center. We have generated 100,000 samples of magnification 
in the Monte-Carlo simulations. We divided the 100,000 samples into 100 datasets of 1000 
samples, and performed the x 2 -test for each dataset, since the same analysis might be possible 
for the observational data of SNAP. Then, we found the significant difference in the rate of 
rejection between n < 0.7 and n > 0.8. This result suggests that we may obtain the 
information about a slope of the density profiles of the central cusp of dark matter using 
MPDF of Type la supernovae. 

A mathematical explanation of the reason why we obtain the gamma distribution is 
unavailable now. In Refs.22), 33), 55), and 25), analytic fitting functions different from ours 
for MPDF are proposed. The relation between our analytic fitting function and those should 
be studied. We note that we need to examine the effect of simplifications carried out in this 
study, before we compare our results with observations. For example, we have assumed that 
there are no spatial correlations of the distribution of the clumps. The effect of the spatial 
correlation of the clumps should be clarified. We leave these issues as our future works. 

The density profile of dark matter halos is a hot topic in astrophysics both theoretically 
and observationally (See e.g. Ref.29)). Our results suggest that, although there are some 
simplifications in comparison with realistic situations of our universe, the statistical gravi- 
tational lensing effect may shed a new light on the observational investigation of the density 
profile of dark matter halos. 
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Appendix A 

Calculation Method 

To calculate magnification factors, we use the multiple lens-plane method. 37 ) We consider 
the "straight" line A from the source to the observer, and we put N lens planes between them 
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so that the straight line A intersects vertically. We denote lens planes from the observer to 
the source sequentially as Si, S2, ■ ■ ■ , ^n- For convenience, we label the source plane that 
is also orthogonal to the line A as Sn+i- On each lens plane, say Si, the lens position 
and the ray position 7^ are specified with respect to the intersection point of A. 

Suppose lens positions (i = 1, . . . , N) are given. When we emanate a ray with a ray 
position 7 X on Si, the lens position on the j-th lens plane Sj (j < N + 1) is recursively 
given as follows 

D j=i 

-rj = ^ri-52 D v*(rYi-Ci), (a-i) 

^ i=l 

where is the DR distance from the i(< j)-th plane to the j-th plane. The right hand 
side of equation fl A- 1 j) also includes ray positions 7 2 , . . . , 7j_i- Thus, we must calculate the 
ray positions 7 2 , . . . , 7j_i to obtain the ray position 7 .. For later convenience, we specify 
the ray positions and lens positions using the following dimensionless quantities: 

Ui : = (A-2) 

Qi ■= jy- (A-3) 

Equation (1A-ljl is written as 

u 3 = u l-'52Pi3 a ( u i-<Ii)> ( A4 ) 

i=l 

where 



DjjD s 
Pij ■ DjD 

and 



.. : = ri2^£ (A- 5) 

3 DjD iS K } 



D s 

OL{ui-q % ) : = — a(7 i -C i )- (A-6) 
Ms 

In the case of the multiple lensing effect, multiple images appear in general. We specify 
the paths from the source to the observer using their ray positions, say v^, on the first lens 
plane. For a path Vm, the magnification is defined as 

= d^v (A ' 7) 

where Am is the Jacobian matrix of the lens mapping: 

du N+1 



A, 



(p) 



dui 

N 



u 1= v 



(p) 



y—^ da(ui - q t ) Out 

^— i OUi OUi 

1=1 



(A- 

Ui=V {p) 
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Observer 



Fig. 12. The definitions of vectors ocj, Cj, and ~fj are given here. 



We obtain the following result. 

Proposition. The magnification probability distribution function for the identical lens objects 
is invariant under the following constant rescaling of the linear extent R and the mass M L 
of the lens object 

R — > CR and M L ^C 2 M L , (A-9) 
if the surface mass densities of the lens objects are given in the form of Eq. [3^4\). 



Proof. The parameter £o becomes C£o through the changes of the parameters (1A-9I) due 
to its definition (12-101) . Hence, we find from Eq. (12- 12[) that the density p{y,z) of lenses in 
(y, z) space is unchanged through the replacement (1A-9I) . if ( is also replaced by C(; hereafter 
we assume so. 

Substituting Eq. (13 -4p into Eq. (12-11) . we have 



a(6 = Const X «j£ / (A.10) 



'12 



For the changes in the parameters (1A-91) . a changes in the manner 

&(£)—> Const x —± I ,V \n 2 — d i 



R 2 J£'<CR l£ ~~ £ 

4M L /• (Z-Ci)F(R/0^ , 



= Ca(C-^), (A-ll) 
where in the second line, we have replaced the variable of integration £' by £ = C^ 1 ^'. 
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The solution of Eq. ( 1A-1I) with the parameters changed as Eq. ( 1A-9I) is denoted by 7^. 
By using Eq. (lA-llj) . the equations satisfied by 7^ are given as 



It is easily confirmed that 7 • = C~fj is the solution for the above equation, where 7 ■ is the 
solution of Eq. AA-IH . Thus, for the changes in the parameters of lens objects (1A-9|) . we have 
the following scaling law 

Uj — >Cuj, (A- 13) 

(see Eq. ( 1A-2I) ). Since the Jacobian matrix A( p ) is invariant for the constant rescaling ( 1A-13I) 
due to definition ( 1A-8I) . we eventually find that the Jacobian matrix A[ p ) does not change 
through the replacement flA-9j) . This result means that the magnification is invariant 
for the change in parameters in the form of (1A-9I) and thus the magnification probability 
distribution function is also invariant for (lA-9h . □ 

The above result guarantees that the MPDF depends on only the parameter R/\/Ml- 
We randomly put lenses so that the distribution of those is consistent with equation (12- 12[) . 
First, we divide the spherical region z < zs into N concentric spherical shells each of which 
is bounded by two spheres z — Zi — Az/2 and z = Zi + Az/2. We take into account 
only the nearest lens in each shell. We find from equation (12- 12[) that in the i-th shell 
Zi — Az/2 < z < Zi + Az/2 there is one point mass on average within the region y < Y i: 
where 

Yi 

p(y,Zi)dy = l. (A-14) 





Therefore we randomly put a point mass within the region y < in the z-th shell. We 
can neglect the lensing effects of the clumps that are so far from the ray that they do not 
affect the magnification. We set an upper bound y max (< Yi) to the lens position, and take 
account of lensing effects due to the lenses in the region y < y max . We have set the value of 
2/max as0 

y max = max ^5, ^ + lj . (AT5) 

Next, we show the method to find paths from the source to the observer: V/ p ). When 
lens positions are given, the dimensionless ray position u N+ i on the source plane Sn+i is 
given as a function of only the ray position U\ on the first lens plane S\\ 

u N+ i = /(ua). (A-16) 



*) We have determined the value of y max as sufficiently large so that < /i >~ D 2 DR / D 2 FL in the case of 
the power law tail lens model. 
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Therefore, to find v^, we solve the equation: 



u N+1 { Ul )=0. (A- 17) 

We find the roots of equation (IA-171) by the Newton- Raphson method. The technical details 
can be seen in Ref.7). 

Appendix B 

Typical Value of the Magnification for the Lens Model (a) 



In the case of the lens model (a), from Eq. (12-11) . the bending angle vector is given 

as 

4Mr 

«(0 = -ffr* (B-l) 
Since the magnification \i due to only one lens plane is given as 

1 Dl 



l> 

we have 



(B-2) 



R 

V = 7 ri (B-3) 

(fit. _ 4M L D L D LS \ 

for the lens model (a) from Eq. (12-21) . Using D$ instead of 2D^ or 2D^s, we have 

H o . (B-4) 

(Ri-M L D S ) 2 
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